Sensitivity of the spherical gravitational wave detector MiniGRAIL operating at 5 K. 
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We present the performances and the strain sensitivity of the first spherical gravitational wave 
detector equipped with a capacitive transducer and read out by a low noise two-stage SQUID 
amplifier and operated at a temperature of 5 K. We characterized the detector performance in 
terms of thermal and electrical noise in the system output sygnal. We measured a peak strain 
sensitivity of 1.5 ■ W~ 20 Hz~ 1/2 at 2942.9 Hz. A strain sensitivity of better than 5 • 10-20#z-l/2 
has been obtained over a bandwidth of 30 Hz. We expect an improvement of more than one order 
of magnitude when the detector will operate at 50 mK. Our results represent the first step towards 
the development of an ultracryogenic omnidirectional detector sensitive to gravitational radiation 
in the 3kHz range. 

PACS numbers: 04.80.Nn, 95.55.Ym,07.07Mp,02.60.Pn 
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I. INTRODUCTION 

The direct observation of gravitational waves (GWs) 
is one of the most challenging tasks for experimental 
physics. After the first detection will be claimed a new 
branch of astronomical observation will begin and grav- 
itational wave observatory will become more and more 
common facilities. A spherical detector is a perfect in- 
strument for an astronomical observatory due to its fea- 
ture of omnidirectionality and polarization sensitivity 
[3, S, H H H . The first ultracryo genic spherical gravita- 
tional wave detectors @, @] are currently completing their 
engineering phase and will soon be operational with an 
expected sensitivity better than 10 _21 77z -1 / 2 at 3kHz. 
We report the results of the first sensitive measurement of 
the spherical gravitational wave detector Minigrail. The 
detector read-out is based on capacitive resonant trans- 
ducers coupled to a superconducting quantum interfer- 
ence device (SQUID) linear amplifiers by means of su- 
perconducting transformers. The two-stage SQUID de- 
scribed in this work is one of the most sensitive amplifier 
ever used on a gravitational wave antenna. It consists 
on a sensor dc SQUID amplified by a Double Relaxation 
Oscillation Squid (DROS). The coupling of the SQUID 
system to the high Q electrical resonators is similar to 
the one developed for the AURIGA detector @ We ob- 
tained a coupled additive energy resolution of 700ft at 
5K in agreement with the expected values calculated 
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from the SQUID parameters using the standard model 
@. The spherical antenna described in this paper is the 
first example of a multimodal resonant detector where 
the five quadrupolar modes of the sphere are read-out by 
three resonant transducers. In this paper we discuss the 
noise contribution and the signal response of one read- 
out channel and estimate the detector sensitivity when 
the detector will operate at 50 mK with a complete read- 
out. This paper is organized as follows. In section |TT] 
we described the experimental apparatus and in particu- 
lar the read-out system. In section section Hill we present 
and discuss the experimental results. In section lTlI Al and 
section HlIBI we analyzed the electrical system, the noise 
spectra and equivalent temperature of the resulting cou- 
pled oscillators. Finally in section HlI CI we describe the 
calibration procedure and estimate the detector strain 
sensitivity. 



II. SYSTEM OVERVIEW 

A. Sphere and mechanical transducers 

MiniGRAIL is a spherical gravitational wave (GW) an- 
tenna currently under development [||. The antenna is 
a massive sphere in CuAl, has a diameter of 68 cm, a 
mass of about 1.3 ton and the GW sensitive spheroidal 
quadrupole modes have frequencies around 2980 Hz at 
4.2 K. The alloy CuAl6% has been chosen because of the 
high quality factor (Q ~ 10 7 at low temperature), high 
sound velocity (Vs — 4100 m/s) and a sufficient thermal 
conductivity, which allows to cool a 1.3 ton antenna to 
a temperature below 100 mK [nj. The ultimate goal 
is to operate MiniGRAIL at a thermodynamic tempera- 
ture of 20 mK, equipped with six transducers coupled to 
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nearly quantum limited double-stage SQUID amplifiers 
[111, [12] ■ The sphere is suspended from the centre with 
a gold-plated copper rod 20 mm in diameter. The rod 
is connected to the last mass of the mechanical vibra- 
tion isolation system which consists of seven mass-spring 
stages suspended with stainless steel cables from three 
absorbers, each consisting of a stack of rubber and alu- 
minium plates. A detailed description of the detector 
mechanics and cryogenics can be found in [l3j 

We used capacitive transducers to read-out the 
spheroidal modes. They consist of a closed membrane 
with a load mass in the centre. The electrode is made 
of a thin CuAl plate placed in front of the resonating 
mass. To obtain a small gap between resonator and elec- 
trode, we proceeded as follows. The resonator and the 
electrode are lapped and polished to get a smooth flat 
surface. Further a clean Kapton foil of a thickness equal 
to the desired gap is placed between the electrode and the 
mass. Finally a small amount of glue is added between 
the electrode and the support springs. A load is applied 
on top of the electrode in order to make a compact as- 
sembly After the drying period of the glue the Kapton 
foil is removed. This technique was shown to be reliable 
and reproducible. Gaps of the order of 20/ito could be 
obtained and voltage bias as large as 500 Volt could be 
applied without discharging. 

In table U we summarise the features of the three trans- 
ducers. Each resonator has an effective resonant mass of 
about 200 g and is tuned mainly to three different spher- 
ical modes. 

Closed membrane transducers 

transducer 1 transducer 2 transducer 3 



mass [Kg] 


0.205 


0.153 


0.150 


Ct [nF] 


1.17 ±0.01 


0.70 ±0.02 


1.20 ±0.01 


gap [ft m] 


20 ± 2 


35 ±4 


25 ±2 


fres [HZ] 


2863 ± 5 


2850 ± 5 


2878 ± 5 


Q at 300 K 


1.0 ■ 10 4 


1.1 ■ 10 4 


1.0 ■ 10 4 


Q at 77 K 


3 • 10* 


4.8 ■ 10 4 


2.0 • 10 4 



TABLE I: Properties of the three CuAl6% closed membrane 
transducers. The transducers resonance frequencies has been 
estimated from the tuning procedure at room temperature as 
described in um 



B. The calibrator 

The calibrator is a capacitive resonant transducer. The 
impedance of the calibrator, biased with an electric field 
E ca i , for each mode of resonance u m , is given by 
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The real part is used to estimate the energy of the mode 
and can be derived as follows 
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on the calibrator method described here particularly in- 
teresting because it is free from systematic errors, at least 
to estimate the temperature of the mode. 



C. The read-out system 

Two transducers, named transducer 1 and transducer 
2 were coupled respectively to the two-stage SQUID and 
to a single stage commercial Quantum Design SQUID. 
The third resonator, transducer 3, was coupled to a room 
temperature FET amplifier and was used for diagnostic 
and calibration purposes. Here we describe the perfor- 
mance of transducer 1 coupled to a two-stage SQUID 
system based on a DROS as an amplifier and a 

Quantum Design dc SQUID as sensor SQUID. [Hj]. The 
complete read-out circuit of transducer 1 is shown in fig- 
ure [1] 

The impedance matching between the transducer ca- 
pacitance and the SQUID input coil is achieved by using a 
high-Q superconducting transformer. The electrical res- 
onance of the transformer is not tuned to the mechanical 
modes, so the impedance matching is not optimal. The 
transformer coils were made of Nb wires and enclosed 
into a double superconducting shield. The bias circuit 
and the decoupling capacitor have been housed into a 
separate compartment of the superconducting shielding 
box. The measured transformer primary and secondary 
coil inductances are respectively L p = 0.3595 ± 0.005 H 
and L s = 2.1 ± 0.2 fxH. The mutual inductance between 
the coils was M = (2.5 ± 0.5) • 10~ 5 H and the coupling 
factor a P: 3 = 0.43. 

We use a decoupling Teflon commercial capacitor d = 
220 nF and SMD bias resistors for a total resistance of 
Rbias = 13 Gil at 4.2 K. The final intrinsic electrical 
quality factor, after connecting the decoupling capacitor 
Cd and the bias resistor Rbias, is equal to Q = (1.8 ± 
0.1) • 10 5 at 4.2K. 

The transducer has been assembled with a gap of about 
20 ± 2 /im and a capacitance C = 1. 17±0. 01 nF measured 
at room temperature. This has been done in order to 
keep the electrical mode separated from the mechanical 
ones. The electrical mode resonance frequency measured 
on the antenna at 5K was 8120 Hz. 

Two-stage SQUID systems are developed in order to 
reduce the noise of dc-SQUID amplifiers, which is nor- 
mally limited by the room temperature electronics (l5| 
H3 UH [3 ■ When used in the transducer chain for reso- 
nant gravitational waves detectors, they can improve the 
detection sensitivity by orders of magnitude Q . The sys- 
tem described here differs from other two-stage SQUIDs 
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dc SQUID 




FIG. 1: Electro-mechanical scheme of a spherical antenna 
with mechanical resonator and capacitive transducer cou- 
pled to a SQUID through a superconducting matching trans- 
former. 



used in GW experiments since it uses a DROS as an am- 
plifier SQUID 13. A DROS has a large flux-to- voltage 
transfer function which allows direct read-out of the sig- 
nal. Direct read-out simplifies multichannel read-out as 
needed in spherical gravitational wave detectors. 

The two-stage SQUID system we developed is based 
on a configuration reported in [l5|, H(| • A standard dc- 
SQUID chip manufactured by Quantum Design (QD) [3l[ 
was chosen as sensor SQUID because of a larger input 
inductance with respect to the dc-SQUID described in 

The dc-SQUID is biased at a constant voltage by 
means of a resistor Rbias = 1-50. The current through 
the sensor SQUID is modulated by an applied signal flux 
Q S ig and is fed through the input coil of the DROS. 
The total addi tive flux noise at 4.2 A' is \fS& = 1.60 ± 
0.02 fi^a/y/Hz with input coil open and \/S$ — 1.10 ± 
0.02 /i<&o/\/Hz with input coil superconductily shorted. 
This corresponds respectively to an intrinsic uncoupled 
energy resolution of e = S<$>/2L sq = 650 ± 15 H and 
e = 320 ± 15 h. This is in agreement with the expected 
values calculated from the SQUID parameters using the 
standard model In order to avoid instability in the 
SQUID-resonator system we implemented a capacitive 



cold damping network in the feedback line. Damping net- 
work has been first investigated by Stevenson [21| , using 
a phase-shifted inductive feedback, and by Vinante Q, 
who made use of a capacitive network. The two-stage 
SQUID coupled to a high quality factor electrical res- 
onator showed the same performances [14| . We estimate 
at 4.2 K a SQUID noise temperature T N = 100 ± 30/jK 
and a noise number N = 730 ± 100. The additive coupled 
energy resolution was 650?i and 320H respectively at 4.2 
and 2.1 K. 



III. EXPERIMENTAL RESULTS 
A. Electrical system and noise spectrum 

The additive noise level of the two-stage SQUID cou- 
pled to the transducer mounted on the sphere, was com- 
parable with the one measured with the SQUID with 



open input. When operating without the cold damp- 
ing network, the minimum wideband flux noise observed 
with the SQUID was of ~ 2.7 ^ /VWz. When the cold 
damping was active, we measured a additive wideband 
flux noise of 1.67 ± 0.03 /i$o/V Hz. It corresponds to 
an additive coupled energy resolution of 730 ± 100 h. As 
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FIG. 2: Flux spectral density measured at the SQUID output 
with transducer bias at 200 Volt. All the expected 8 modes 
of the system are visible in the spectrum. 

shown in figure [H all the expected 8 modes of the system 
are visible in the spectrum. To establish which modes 
were strongly coupled to the transducer with the two- 
stage SQUID system, we proceeded as described in the 
next section. 



B. Equivalent temperature of the coupled 
oscillators 



The best and most direct way to estimate the noise 
of the system is to measure the input impedance of the 
transducer coupled to the read-out SQUID amplifier as 
described in [22|, [H[. However, in this experimental test 
we did not implement in the matching transformer the 
necessary calibration coil. To estimate the temperature 
of the modes we can proceed as follows. The power spec- 
tral density at the output of the calibrator, when the 
back-action contribution of the room temperature am- 
plifier is negligible, is 



Sv,cal = ik B T eq Re(Z(uj)), 



(3) 



where Re{Z{u>)) was derived in Eqs. ([T]), T eq is the equiv- 
alent temperature of the modes and we are considering 
monolatcral spectra. At low temperature this value is too 
small to be measured with the room temperature FET 
amplifier. However, by exciting each modes at resonance 
with an auxiliary piezo-electric transducer (PZT), we can 
increase the signal at a level that can be measured by the 
FET amplifier. Then we read both the responses of the 
FET amplifier and the two-stage SQUID coupled to the 
resonator biased at a voltage Vb. sq - We assume that the 



4 



voltage Vcal (u>) at the output of the FET is proportional 
to the voltage Vsq(u>) at the output of the SQUID, i.e. 
Vcal = A ca i ySq VsQ((jj)- This is true at resonance. If we 
now measure the power spectrum Sy,SQ at the SQUID 
output, when the modes are not excited, we can evaluate 
the equivalent temperature of the modes using Eqs. (|3j). 
where we substitute 



~>V,CAL 



A. 



,S 



V.SQ- 



(4) 



Here we assumed that the system response is linear in 
the whole range. Linearity has been checked for different 
excitation voltage. We found a linear behaviour within 
10%. 

To estimate the equivalent temperature of the modes 
we measured the variance a 2 of the stochastic process 
with a spectral noise Sy,SQ at the output of the two-stage 
SQUID. From Eqs. Eqs. © and Eqs. ©, defining 



m m Q m 



the variance can be written as follows 



2tt A 2 



(w m +A) 



dtu 



l, sq J(u m -A) {ujI - u 2 ) 2 + 



(5) 



from where we obtain the relation which links the equiv- 
alent temperature of the mode to the variance of the 
stochastic process with spectral noise given by Sy,SQ- 
We have 



T, 



eq 



(/ ,3 A2 

^ra^cal.sq ^2 

2kBA m flQ m m 



(0) 



where A 2 al sg comes from the calibration as described in 
this section, and ^4 m ,o is estimated from the Lorentzian 
curve fitting of each resonance of the real part of the 
calibrator impedance. 

To estimate a m a lock-in amplifier is used with the 
reference set at the resonance frequency of the mode. 
The lock-in amplifier output magnitude r and angular 
phase <p are then sampled at regular time intervals. The 
amplitude decay time constant of the lock-in amplifier is 
chosen equal to the sampling time t s = Tik- To observe 
the free evolution of the mode m the lock-in amplifier 
time constant is chosen smaller than the time constant 
of the mode, Tik < T m , but large enough that the lock- 
in amplifier works as a bandpass filter and makes the 
contribution of the broadband noise of the SQUID and 
the tails of the neighboring modes negligible. The mean 
square amplitude (r 2 ) of the lock- in amplifier input signal 
magnitude, equal to the variance of the total narrow- 



band noise V 2 b , is given by [24 1 



Swb \ 

2~Vk~J ' 



(7) 



where Cq is the variance of the power spectral density out- 
put, and S w b is the power spectral density of the SQUID 
wideband noise. Generally the factor T TO /(r m + Tik) ~ 1 
in our case. 



The stochastic process r 2 is the sum of two indepen- 
dent Gaussian processes, the in-phase and quadrature 
lock- in amplifier output. If the signal is absent or in gen- 
eral if its average contribution is negligible with respect 
to the noise, the variable r 2 will have the exponential 

2 

distribution F(r 2 ) = r^—e - ^. 

The estimate of Oq , is then performed by sampling the 
magnitude r at regular time intervals, with sampling time 
t s << r m . A subset of data is created by extracting 
a data point every resonator time constant r m in order 
to get uncorrelated samples. After a large number of 
samples is collected, a histogram N(r 2 ) is built, where 
N is the number of samples in a given bin around r 2 . 
The histogram is fit with the exponential distribution 
described above and the mean square amplitude (r 2 ) is 
then extracted as fitting parameter. 

In the absence of excess or amplifier back-action noise, 



the quantity 



is proportional to the thermal vibra- 



tional energy in the antenna mode. The constant of pro- 
portionality a m was used to rescale the recorded values of 
(r 2 ) to antenna energy. For the two most coupled modes 
at frequencies 2943 Hz and 2985 Hz we found the cali- 
bration factor a m to be a 2 943 = (7.0 ± 1.5) • 10 8 [K/V 2 ] 
and a 2985 = (1.1 ± 0.2) • 10 9 [K/V 2 ] respectively. 

Graphs (a) and (b) in figure[3]show the energy distribu- 
tion estimated for the modes at frequencies 2943 Hz and 
2985 Hz during three hours of acquisition. The equiva- 
lent temperature for both modes is obtained by fitting 
the exponential distribution. 

The slope of the distribution corresponds to a temper- 
ature of 7.0±2K for the mode at 2943 Hz and 9±2 K for 
the mode at 2985 Hz, the error arising mainly from the 
calibration uncertainty. The equivalent temperatures are 
consistent, within two sigma, with the thermodynamic 
temperature of the sphere. No significant difference is 
observed in the equivalent temperature of the mode be- 
tween night and day acquisitions. 



C. Force calibration and strain sensitivity 

When the mode m is excited at resonance uj m with 
an energy given by ^ksT m , the power spectral density 
measured at the SQUID output is 



S v 



Ak B T m Re{Z(io m )) 
sq A 2 

cal,sq 



V 2 



Hz 



(8) 



where Re(Z(ui)) is the real part of the calibrator 
impedance. 

The force power spectral density of a mode excited at 
a temperature T m , when back action is negligible like in 
our case, is given by 



S 



FF,m — 



4kBT m m e ff,, 



TV 2 

TTz 



(9) 
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FIG. 3: Exponential distributions of the mean square amplitude (r 2 ) for the modes at 2943 Hz and 2985 Hz. The variance 
obtained from the fit of the exponential distribution gives an equivalent temperature of the two modes of 7.0 ± 2 K and 9 ± 2 K 
respectively 



where m e ff >m is the effective mass of the mode m, which 
can be estimated from the tuning curves of the sensor 
transducer if the bias voltage and the gap are known. 
[3^ From the ratio of Eqs. © and Eqs. © we get 



sq 



' l eff,m^mA cal ^ q 
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From the calibrator impedance measurement, we have 



Re{Z{u>)) 



An 



(11) 



where w mjC and Q m ,c were measured during the calibra- 
tion process and were different from the resonance mea- 
sured during the acquisition of the mode spectra because, 
in the later case, the calibrator was not charged. A m c is 
obtained from the Lorentzian fit of each mode of the real 
part of the calibrator impedance. 

At each mode resonance Eqs. (|TU|) becomes 



sq 



>FF,i 



) 2 

■> m.c 



XL 
N 2 



(12) 

We now turn to analyze the detector transfer func- 
tions. The relations derived so far are valid at resonance. 
As a first approximation, one can consider the transfer 
function Gsq,f as a product of poles and zeros where the 
poles are derived from the polynomial fit of the SQUID 
noise spectrum and the zeros are chosen to fit the mea- 
sured amplitude at resonance given by Eqs. (|12D . 

The transfer functions for an applied calibration signal 
becomes 



(13) 



In the equation above N p > N r and H m ^ ca i(uj) is a force 
calibration constant which has been experimentally de- 
termined from the calibration measurement at resonance. 

All the terms included in Eqs. (fl"2|) and Eqs. (Tl~3"|) are 
derived experimentally from the calibration, from the 
tuning curves and from direct measurement of the modes 
quality factor. The transfer function is experimentally 
measured at resonance and only approximated out of res- 
onance. 

We remark again that quality factors and resonance fre- 
quencies are different when measured during calibration 
(w TOlC , Q m ,c) and during acquisition of the noise spectra 
(Lu m ,s,Qm,s), due to the bias voltage in the calibrator. 
This complication arises from the fact that the calibra- 
tor was also coupled to the quadrupolar modes of the 
sphere. This effect is included in the estimate transfer 
function of Eqs. JT3]). 

A small parenthesis needs to be opened here. While 
for a bar detector it is relatively straight forward to re- 
late the strain produced by a calibrator located on one 
of the bar faces with the strain from a gravitational wave 
signal, the same cannot be said for a spherical detector. 
A calibration hammer stroke excites a linear combination 
of the five spheroidal modes depending on the position 
on the sphere surface. When only one calibrator is used, 
like for the Minigrail test run described in this paper, one 
only calibrates the detector for a particular set of forces 
F m applied to each spheroidal mode. Such a combination 
of forces might not always represent a GW excitation. In 
order to fully calibrate the detector one needs to gener- 
ate calibration forces from a set of 5 or more calibrators 
located at different positions on the sphere surfaces. A 
detailed procedure to calibrate a spherical detector is de- 
scribed in [23l | . 

In the experiment described here, the calibration signal 
is generated by a piezo located at position (9 = 18°, <f> = 
135°). By using Eq. 45 in [23|, derived previously in 
[H, d, [25| , one finds that a piezo in such a location excites 
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a combination of the 5 spheroidal modes given by the vec- 
tor (0,0.13,-0.49,0.5,-1), normalised to the maximum 
value of its elements. One can see that such a combina- 
tion is equal, within 20% tolerance, to the one generated 
by a circularly polarised gravitational wave coming from 
direction (9 = 20°, cj) = 135°). Following [23| one can 
see that such a wave direction is almost optimal for a 
transducer configuration described in this paper. 

MiniGRAIL strain sensitivity with three transducers 
coupled to the spheroidal modes, but only one used 
for the read-out, is shown in figure 0] The read-out 
transducer was biased with an constant electric field of 
E = 10 7 Volt /m and the sphere thermodynamic temper- 
ature was 5.2 K. As discussed above, the experimental 
strain curve gives and estimation of the detector sensitiv- 
ity only for a particular combination of spheroidal modes 
corresponding to a gravitational wave coming from the 
(9 = 20°, = 135°) direction. 




frequency [Hz] 



FIG. 4: MiniGRAIL strain sensitivity at 5 K with three 
transducer placed on the sphere, but single transducer read- 
out. The sensitivity has been estimated for a particular com- 
bination of the 5 spheroidal modes as derived by exiting 
the mode using a piezo transducer (PZT) located at (9 = 
18°, (f> = 135°). The dashed gray curve shows the strain sen- 
sitivity calculated for a simulated hammer stroke excitation 
from the PZT location using the electro-mechanical model 
described in [5^]. For the simulation we used the detector 
parameters discussed here in the text. A better matching be- 
tween the experimental data and the simulation is obtained 
when the simulated hammer stroke is given at the position 
(9 — 27°, <f) = 135°) (continuous gray curve). The other 
curves shows the contribution of the thermal noise (dashed 
dark line) and the back-action noise (dotted line) to the strain 
sensitivity 

We obtained a peak strain sensitivity of (1.5 ± 0.6) • 
IQ-^Hz- 1 ' 2 at 2942.9 Hz and a strain sensitivity of 
about 5 • IQ-^Hz- 1 / 2 over a bandwidth of 30 Hz. This 
corresponds to a strain amplitude of h ~ 2.5 • 10~ 18 
at 3 kHz for a burst signal of 1 ms [26| . For a sphere 
of 68 cm in diameter like the one of MiniGRAIL, it is 
equivalent to a displacement sensitivity, at 3 kHz, of 



1.6 ■ 10 _19 m. When optimal filtering is applied to the 
output signal, the detector is sensitive to burst signals 
with an impulse energy of about Tjv ~ 50 mA'' as can 
be derived calculating the noise temperature using the 
experimental data. This sensitivity would be enough to 
detect supernova explosions in our galaxy. The calcu- 
lated sensitivity curves in figure 01 arc obtained using 
the model described in [23| . The dashed gray curve shows 
the strain sensitivity calculated for a simulated hammer- 
stroke excitation applied to the same point on the sphere 
where the PZT is placed, i.e. (9 = 18°, = 135°). The 
simulated signal describes reasonably well the strain sen- 
sitivity 

The best fitting strain sensitivity has been obtained for 
a simulated hammer stroke applied at the sphere surface 
point (9 = 27°, = 135°). The result is shown with the 
continuous gray line. In this case the third mode is more 
excited. The agreement with the experimental data is 
rather impressive considering the amount of fitting pa- 
rameters involved in the simulation. The difference of 
about 7° in the angle 9 of the experimental and simu- 
lated excitation position could be explained considering 
the fact that a spherical detector with only three trans- 
ducers in the position 1,2 and 5 like the one considered 
here is far from being symmetric. In [2^ | the authors had 
to perform a rotation of spheroidal mode reference frame 
as well to be able to explain their experimental results. 
In their case the transducers were not as massive as here 
and, above all, they used six transducers positioned in the 
symmetric TI configuration. The resulting mixing of the 
spheroidal modes could explain the discrepancy between 
the measured and the simulated sensitivity curves. In or- 
der to address more accurately this issue one should place 
on the sphere at least 5 calibrators to measure the detec- 
tor response to each of the 5 spheroidal modes. When 
six transducer are fully operating with comparable sensi- 
tivity, a single calibrator is enough to fully calibrate the 
detector as shown in [H, [27j . 

The contribution to the strain sensitivity of the well 
known noise sources are plotted in figure[Has well. At res- 
onance, the sensitivity is limited by mechanical thermal 
noise of the transducer mass. Out of resonance the sen- 
sitivity is limited by the SQUID additive current noise. 
The back action noise of the SQUID is about an order of 
magnitude smaller. The electrical thermal noise of the 
superconducting transformer, not shown in the graph, is 
negligible because the electrical mode is well decoupled 
from the mechanical ones. 

From the measurement of the variance of the most cou- 
pled modes and the simulated data we can conclude that, 
within the experimental accuracy, the Minigrail peak sen- 
sitivity is currently limited by the thermal noise gener- 
ated by the transducer mass. Some of the modes, how- 
ever, show excess noise whose origin is difficult to ad- 
dress. A better estimate of the transfer function of each 
spheroidal mode is necessary in order to fully characterise 
the detector. The measurements presented here have to 
be considered as a first test bed for the following engi- 
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neering and science runs of Minigrail. The sensitivity is 
expected to improve of at least one order of magnitude 
when the detector will operate at 20 mK. 

In figure [5] the measured strain sensitivity is shown 
together with predicted sensitivity for possible future de- 
tector configurations. A polynomial fit of the strain sensi- 
tivity is shown as well. The poles and zeros obtained from 
the fit can be used to build the matched filters for GW 
detection [23[. The expected strain sensitivity of Mini- 
grail is shown for the detector operating at T — 50 mK 
with the same three transducers configuration presented 
here. The figure shows the expected sensitivity for a Min- 
igrail II. In this configuration Minigrail operates with 
6 capacitive transducers placed in the TI configuration 
where the electrical modes are coupled to the mechanical 
ones. We consider T/Q = 2.5 x 10~ 8 A" and a SQUID 
coupled energy resolution E coup i e d = 70h. Those val- 
ues are achievable within the current technology. See 
[23l | for a recent review. The sensitivity of Minigrail op- 
erating at the quantum limit is also shown. Minigrail 
can reach a peak sensitivity of about 6 • 10~ 2i Hz~ 1 / 2 
and a bandwidth larger than 400 Hz at a sensitivity of 
1 ■ IQ-^Hz- 1 / 2 . 



IV. CONCLUSIONS 



We have operated at 5 K a spherical resonant de- 
tector equipped with a capacitive resonant transduc- 
ers coupled to a two-stage SQUID amplifier. Our two- 
stage SQUID amplifier is one of the most sensitive am- 
plifier employed so far on a GWs resonant detector. 
We measured an additive coupled energy resolution of 
700 ± 100 h at 5K We reach a peak strain sensitivity of 
1.5 ■ lO" 20 ^" 1 / 2 at 2942.9 Hz. A strain sensitivity of 
better than 5 ■ 10~ 20 -ff2i -1 / 2 has been obtained over a 
bandwidth of 30 Hz. We expect an improvement of more 
than one order of magnitude when the detector will oper- 
ate at 50mK. This result should be considered as the first 
step towards the realization of an ultracryogenic spheri- 
cal gravitational wave detector. 
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FIG. 5: The measured strain sensitivity of MiniGRAIL is 
shown together with predicted sensitivity for future detector 
configurations. The continuous gray line is a polynomial fit 
of the measured strain sensitivity. The dashed line shows 
the sensitivity for the detector operating at T = 50 mK 
with the same three transducers configuration presented in 
this paper. The dot-dot-dashed line (Minigrail II) shows 
the sensitivity achievable with available technology , namely 
T/Q ~ 2.5 x 10 -8 K and SQUID energy resolution E = 70/L 
The dot-dashed curve gives the sensitivity for a quantum lim- 
ited detector (Minigrail QL) with T/Q ~ 1 x 10~ 9 K . 
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